Non monotonic velocity dependence of atomic friction 
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We study the velocity dependence of the frictional force of the tip of an atomic force microscope as 
it is dragged across a surface, taking into account memory effects and thermal fluctuations. Memory 
effects are described by a coupling of the tip to low frequency excitation modes of the surface in 
addition to the coupling to the periodic corrugation potential. We find that when the excitation 
mode frequency is comparable to the characteristic frequency corresponding to the motion of the 
tip across the surface, the velocity dependence of the frictional force is non monotonic, displaying a 
velocity range where the frictional force can decrease with increasing velocity. These results provide 
theoretical support for the interpretation of recent experiments which find a frictional force that 
decreases with velocity on surfaces covered with a monolayer. 
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I. INTRODUCTION 

In an atomic force microscope, the motion of the tip as 
it is dragged across a substrate provides an efficient way 
to probe atomic friction of surfaces. In fact, the moving 
tip can be regarded as a single asperity which determines 
the frictional force between macroscopic surfaces 1 . Under 
a constant load, the low-velocity motion of the tip on a 
surface exhibits stick-slip behavior as the moving tip hops 
over the corrugation potential defined by the substrate. 
For a quantitative description of friction at the micro- 
scopic level, theoretical understanding of the behavior of 
the tip on the substrate under different conditions is re- 
quired. Recent studies have argued that due to thermal 
fluctuations, the lateral frictional force increases logarith- 
mically with increasing sliding velocity^— . However, it is 
also found experimentally that the frictional force can 
have an opposite behavior, where it decreases with the 
sliding velocity depending on the ambient conditions 5 
and the nature of the substrate^. The origin of this 
decreasing velocity dependence is not fully understood 
and it is possible that distinct effects need to be taken 
into account to explain different experimental conditions. 
Current explanations of this behavior for point contact 
friction are based on extensions of the simplest effec- 
tive model, the generalized Tomlinson model 3,8 ' 9 , taking 
into account additional effects such as dissipation due 
to deformation of the AFM tip£ or the inclusion of addi- 
tional time scales for reorganization of the surface during 
sliding 5 . In particular, the interpretation of recent exper- 
iments on surfaces covered with a monolayer in terms of 
an additional time scale due to surface restructuring^ still 
lacks detailed theoretical support. 

In this work, we study the velocity dependence of the 
frictional force within a simple theoretical model taking 
into account the coupling of the AFM tip to excitation 
modes of the surface in addition to the coupling to the pe- 



riodic corrugation potential. We model the generic exci- 
tation modes by damped harmonic oscillators analogous 
to the model of an adsorbate coupled to phonon modes of 
the substrate^. We find that when the excitation mode 
frequency is comparable to the characteristic frequency 
of the particle moving across the surface potential, the 
velocity dependence of the frictional force is non mono- 
tonic displaying a velocity range where the frictional force 
can decrease with increasing velocity. These results pro- 
vide theoretical support within a simpler model for the 
interpretation of recent experiments^ which find a fric- 
tional force that decreases with velocity on surfaces cov- 
ered with a monolayer. 



II. MODEL 

The lateral motion of the tip can be described by a 
model of a single particle in an external one-dimensional 
potential representing the substrate and coupled elasti- 
cally to the moving support. The equation of motion for 
the particle is given by 3 

dV(R,x) . ,. . 

mx = ; mrfx + j(t) (1) 

ax 

where the noise variables f(t) satisfy 

< /(*)/(*') >= 2 mi k B T5{t - t') (2) 

and V(R, x) — V p (x) + V s (x,R) is the total potential 
including the interaction of the tip to the substrate with 
periodic potential 



V p (x) =V [1-cos{2tt-)} 



(3) 



and the clastic interaction between the tip and the mov- 
ing support 

1 



V s {x,R) = -k(x-R(t)Y 



(4) 
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where R(t) is the position of the support. 

The friction damping parameter 7 in Eq. (JT|) results 
from the coupling of the tip to the surface excitations 
with a much shorter time scale than the motion of the tip 
(e.g. electronic excitations) which results in a (5-function 
correlated random force f(t) with instantaneous Marko- 
vian damping without memory as described by Eq. 
The microscope support is moving at constant velocity, 
R(t) = vt, where v is the support velocity. The average 
force on the support due to the tip motion is then given 
by 

F = k<R(t)~x> (5) 

This is regarded as the frictional force with the substrate 
since without damping we should have F = when av- 
eraging out the negative and positive force region of the 
periodic potential. For sufficiently large frictional damp- 
ing stick and slip motion is expected. As shown in Ref. 
H, using this model the frictional force F increases with 
increasing velocity v. At small v, it increases logarithmi- 
cally but ultimately one would expect F to be just jv, 
for sufficiently large v. 

We extend the model described by Eq. (TTJ) by allowing 
the tip to couple to other excitations in the substrate 
which have a finite time scale. These other excitations are 
modeled by damped simple harmonic oscillators as done 
in Ref. [l^ • These excitations can be labeled by the mode 
index A = {i, I) where i = 1, 2 labels the symmetry group 
and I runs from 1 to N/2 where N is the total number 
of modes. Each mode is characterized by a displacement 
variable u\, which couples to the tip motion. The total 
Hamiltonian can be written as 

A 

+ D^ + ^K + ^M) 2 ] (6) 

The third term in the above Hamiltonian represents 
a counterbalance term. This is added so that when the 
tip moves very slowly relative to all the time scales given 
by 2tt/lo\, then at any instantaneous position of x, u\ = 
—W\ and the actual effective potential felt by the tip 

is V p (x) — ^2 2 X ^A w hich is the adiabatic potential. 
However, when the tip is moving very fast relative to such 
time scales, then the substrate has no time to respond, 
and one recovers the rigid substrate potential V p (x). The 
explicit form of W\(x) is chosen as 

2iraV 

Wi 1 = = sinyiirx/ a) 

W 21 = 2 ™ V ° cos(2irx/a) (7) 

The coupling a in Eq. ([7]) depends on the frequency oj of 
the mode. This is chosen here to have the form a(u>) oc uj 



so that the total effect of a finite number of damped oscil- 
lators corresponds to a continum distribution^ of excita- 
tion modes with a density of states p(uj) oc uj 2 , analogous 
to that for 3-D phonons in the Debye model. The equa- 
tions of motion for the tip and the mode displacement 
variables can be written in a dimensionless form as 

it = —2irsin(2irx) + k(R(t) — x) — 2na cu\f\(x)u\ 
-7± + /(t) 

«A = ojfu x - a a uJxfjx(x) - 7a"a + r\(t) (8) 

where the additional noise variables r\ for the damped 
harmonic oscillator also satisfy a fluctuation-dissipation 
relation 

< r A (i), r A (i') >= 2 lx TS{t - t')6 x ,X' (9) 

Without loss of generality, we can set m = M for simplic- 
ity. In Eq. ©, /i,/ = cos{2ux), gij = sin{2-Kx), f 2 ,i = 
sin(2'Kx), <?2.z = cos(2nx ) and a a = 2nV /(a v / N/2). 
The lattice constant a is taken as the unit of length, V D 
the unit of energy, and \J ma 2 /V the unit of time. 

III. NUMERICAL RESULTS AND DISCUSSION 

We have performed numerical simulations of the cou- 
pled dynamical equations ([8]) using Brownian molecular 
dynamics^. The parameters for the spring and micro- 
scopic damping were fixed to k = 3 and 7 = 6 and for 
damped harmonic oscillators, 7^ = 2w\ and a = 2.5. 
Time steps ranged from dt = 0.001 to 0.005. Numer- 
ical results for the velocity dependence of the friction 
force obtained from Eq. ([5]) assuming different angular 
frequencies lo for coupling to one odd or even mode are 
shown in Figs 1 and 2. They are qualitative similar and 
show that when the frequency of the excitation mode is 
sufficiently small the velocity dependence of the friction 
force is non monotonic. While at sufficiently small or 
large velocities it increases with velocity there is an in- 
termediate range of velocities where the frictional force 
decreases as the velocity increases. The approximate lin- 
ear decrease in the semilog plot of Figs 1 and 2, suggest 
that this decrease is logarithmic. 

To understand this non monotonic behavior it should 
be noted that there are several important time scales in 
the problem. First, there is a velocity dependent time 
scale, t v = a/v, which is the average time for the tip to 
traverse one lattice period of pinning potential. The pin- 
ning potential introduces a time scale, to = \J ma 2 /V , 
the vibrational period of the tip in the well. In the stick 
and slip regime where the tip is stuck in the well and 
perform many vibrational oscillations before hopping to 
the next minima, we have t v » tg. The time for cross- 
ing the barrier is relatively fast, on the order of to. With 
the coupling to a new excitation mode of frequency a;;, 
another time scale t t — 2u/ui is introduced. We now 
consider what is the effect of varying t v — a/v on the 
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FIG. 1: Velocity dependence of the frictional force of the tip 
moving on a substrate coupled to an excitation mode with 
different angular frequencies u>. Temperature kT/V = 0.3. 
Results are for the odd-mode model. 



FIG. 3: Velocity dependence of the frictional force of the 
tip moving on a substrate coupled to an excitation mode of 
angular frequency to — 0.1 at different temperatures T/Vo- 
Results are for the odd-mode model. 
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FIG. 2: Velocity dependence of the frictional force of the 
tip moving on a substrate coupled to excitation mode with 
different angular frequencies uj. Temperature kT/V = 0.3. 
Results are for the even-mode model. 



frictional force F. In the absence of the coupling to the 
excitation mode, the logarithmically increase of the fric- 
tion force with velocity^ was obtained for velocities up 
to t v ps to. Now consider what happen if we choose an 
excitation mode such that U is several times larger than 
to (uji << loo ) and we vary t v . As far as the intrin- 
sic frictional damping is concerned, when t v approaches 
U, the damping should increase and also becomes non- 
Markovian. Eventually, when t v becomes very large, the 
damping should saturate to a Markovian value. So based 
on this consideration, the effective 7 should increase and 
saturate as a function of decreasing v. The frictional 
force F should exhibit similar behavior. In addition, 



the even mode can lead to a velocity dependent sub- 
strate relaxation that changes the effective potential as 
discussed previously. However, for the present choice of 
parameters, the velocity dependence of the effective non- 
adiabatic frictional force due to the coupling to the new 
excitation mode dominates as evidenced by the results in 
Fig. 1 and 2, showing similar results for coupling to an 
even or odd modes. 

The effect of varying the temperature at fixed mode 
frequency u> is also of interest. As can been from Fig. 
3, the qualitative non monotonic behavior remains but 
there is an increase in the magnitude of the negative slope 
of the velocity-force curve as the temperature decreases. 

We have also performed additional calculations to 
check the effect of including more excitation modes to 
model a more realistic continuum model of substrate exci- 
tations. 6 independent damped harmonic oscillators with 
frequencies in a range [ui min , LJ max ] with w mm = uj max /6 
were included in these calculations. As can been seen 
from Fig. 4, the behavior for the velocity dependence for 
different maximum frequencies uj max is qualitative simi- 
lar to the behavior in Fig. 1 at the same temperature. 
The main effect of the additional modes below Lo max is 
a tendency to saturation of the frictional force at lower 
frequencies. 



IV. CONCLUSION 

We study the velocity dependence of the frictional force 
of the tip of an atomic force microscope as it is dragged 
across a surface using a single particle model, taking into 
account memory effects and thermal fluctuations. In our 
model, there is an additional coupling of the particle to 
excitation modes of the surface modeled by damped har- 
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FIG. 4: Velocity dependence of the frictional force of the tip 
moving on a substrate with excitation modes ranging from 
LOmin to u) m ax- 6 modes equally spaced with Lj min = uj ma x/6 
are used. Temperature kT/V = 0.3. Results are for the 
odd- mode model.. 

monic oscillators^. We find that when the excitation 
mode frequency is comparable to the characteristic fre- 
quency of the motion of the tip across the surface po- 
tential, the velocity dependence of the frictional force is 
non monotonic displaying an intermediate velocity range 
where the frictional force can decrease with increasing 
velocity. Therefore, the sliding behavior on the surface 
depends on the nature the surface through the character- 
istic frequency of such excitation modes. Recently, non- 
monotonic velocity dependence of friction was observed 



experimentally^ on surfaces covered with a monolayer. 
The experimental results were interpreted as resulting 
from the chemical nature of the surface, with friction 
increasing or decreasing with velocity depending on the 
presence of cross-linked H-bonds. It was argued that for 
higher sliding velocity there is not enough time for re- 
ordering to occur and the friction force should decrease 
logarithmically with velocity, while at lower velocities the 
breaking and reordering of the H-bonds contribute a new 
source of friction that increases with decreasing veloc- 
ity. Our results provide theoretical support for this idea 
through an explicit calculation within a simple model. 
The breaking of glassy domains of H-bonds at a criti- 
cal stress and reordering introduces a new time scale, 
which corresponds in our model to the inverse frequency 
of the additional excitation mode of the surface. The 
non-monotonic velocity dependence of the friction then 
follows naturally when the time scale of the motion of 
the tip across the surface becomes comparable to this 
new time scale for the excitation mode. However, since 
we find a non monotonic behavior without invoking any 
glassy behavior, our results suggest that similar behavior 
should be found even on surfaces without disorder when 
additional excitation modes are present. 
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